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Abstract. - The geo and solar magnetic fields have long been thought to be very different objects 
both in terms of spatial structure and temporal behavior. The recently discovered field structure of 
a fully convective star is more reminiscent of planetary magnetic fields than the Sun's magnetic field 
[1], despite the fact that the physical and chemical properties of these objects clearly differ. This 
observation suggests that a simple controlling parameter could be responsible for these different 
behaviors. We report here the results of three-dimensional simulations which show that varying 
the aspect ratio of the active dynamo region can yield sharp transition from Earth-like steady 
dynamos to Sun-like dynamo waves. 
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Introduction. - Observations of the magnetic fields due 
to dynamo activity appear to fall into two categories: 
fields dominated by large-scale dipoles (such as the Earth 
and a fully convective star), and fields whith smaller-scale 
and non-axisymmctric structures (such as the Sun). 
Moreover two kinds of different temporal behaviour have 
been identified so far: very irregular polarity reversals (as 
in the Earth), and quasi-periodic reversals (as in the Sun). 
Since the Earth and the Sun provide the largest database 
of magnetic field observations, these objects have been 
well studied and described in terms of alternative physical 
mechanisms: the geodynamo involves a steady branch of 
the dynamo equations, perturbed by strong fluctuations 
that can trigger polarity reversals, whereas the solar 
dynamo takes the form of a propagating dynamo wave. 
The signature of this wave at the Sun's surface yields the 
well-known butterfly-diagram (Sunspots preferentially 
emerge at a latitude that is decreasing with time during 
the solar cycle). 

Modelling. - Because of their very different natures (liq- 
uid metal in one case, plasma in the other), planetary and 
stellar magnetic fields arc studied by different communi- 
ties. Non-dimensional numbers controling the dynamics of 
the Earth and the Sun, for example, do significantly differ 
(see [2,3]). As a practical matter however, the techniques 
as well as the typical parameters used in numerical stud- 
ies of these two systems are surprisingly similar. To some 
extent this is due to the restricted parameter space avail- 



able to present day computations. The parameter regime 
numerically accessible is rather remote from the actual ob- 
jects. For planetary dynamos the main discrepancy relies 
in the rapid rotation in the momentum equation (charac- 
terized by the Ekman number) , whilst for stellar dynamos 
it relies in solving the induction equation with weak resis- 
tive effects (characterized by high values of the magnetic 
Reynolds number). Yet within this restricted domain, the 
sharply different key characters to both geo [4] and so- 
lar [5,6] magnetic fields have been reproduced. This leads 
us to argue that the important parameter controlling the 
magnetic field behaviour is the aspect ratio of the dynamo 
region (i.e. the radius ratio of the inner bounding sphere 
to the outer bounding sphere). Indeed, in the Earth, the 
inert solid inner core extends to less than 40% of the core 
radius, whereas in the Sun, the radiative zone fills 70% 
of the solar radius. One expects the convective zones of 
stars and planets to have all possible intermediate aspect 
ratios, even extending to fully convective spheres. 

In order to isolate and understand this purely geometri- 
cal effect, we have carried out three-dimensional numerical 
simulations of self-excited convective dynamos in which 
the domain aspect ratio was slowly varied, with all other 
parameters held constant. The governing equations as well 
as parameter regimes used here were originally introduced 
for a geodynamo reference calculation [7]. The only dis- 
tinction being the use of stress-free boundary conditions 
on the outer sphere of the domain, while imposing no-slip 
boundary conditions at the bottom of the convective re- 
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Fig. 1: Time evolution of the radial magnetic field averaged in longitude (for an aspect ratio of 0.65). The initial dipole field 
survives for a few diffusion times, and then vanishes to yield a butterfly-like diagram. 



gion. This choice was made in order to create a strong 
shear at the base of the model, and thus try to mimic the 
solar tachoclinc [8]. The inner sphere is here assumed to 
be insulating, and we use differential heating. The gov- 
erning equations are in non-dimensional form: 

E [d t u + (u • V)u] = - Vtt + EAu-2e 2 xu 

+ RarO + Pm- 1 (V x B) x B,(l) 



d t B = V x (u x B) + Pm" 1 AB , 
d t 6 + (u • V)(0 + T t ) = Pr" 1 A6 , 
Vu = V-B = 0, 



where 
E = 
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All simulations reported here were performed keeping the 
following parameters constant E = 10~ 3 , Ra = 100, Pr = 
1, Pm = 5 . The above system is integrated in three- 
dimensions of space (3D) using the Parody code [9] . 

When the inner (non dynamo generating) body oc- 
cupies less than about 60% of the convective body in 
radius, the flow generates a dipolar field, very similar to 
that of the Earth. It features patches of intense flux at 
high latitudes and some reversed patches at low latitude, 
similar to the ones revealed by a downward continuation 
of the Earth's field to the Core-Mantle boundary [10]. 
This strongly dipolar solution becomes unstable with a 
further increase of the aspect ratio. For an aspect ratio of 
0.65 -close to that of the Sun- the strong dipole is first 
maintained and then strongly weakens, but dynamo ac- 
tion continues in a different form: that of a wavy solution 
with quasi-periodic reversals (Fig. 1), reminiscent of some 
aspects of the solar magnetic field behavior. Drifting 
features can be observed both on the radial field at the 
surface of the model (Fig. 1 & 2b) and on the azimuthal 
(east- west) field below the surface of the model (Fig. 2c) . 
Due to the complex nature of these fully tri-dimensional 
simulations, many waves can co-exist. Some of the 
dominant structures appear to propagate toward the 



equator; others propagate poleward. Reversed waves are 
also observed at the surface of the Sun at higher lati- 
tudes [11]. Let us stress however that the model cannot 
be expected to capture all the features either of the geo 
or solar magnetic fields. In particular due to the param- 
eters regime and the lack of stratification in our modelling. 

Physical interpretation. - In order to investigate the 
physical mechanisms associated with these waves, we have 
performed some kinematic simulations. During the course 
of the simulation the Lorentz force was suppressed. The 
wavy nature of the dynamo field was unaltered by this 
modification. This rules out the possibility of an inter- 
pretation in terms of pure Alfven waves or Alfven waves 
modified by rotation (so called MC or MAC waves) , which 
both require the back-reaction of the Lorentz force. Of 
course, suppressing the Lorentz force is not without con- 
sequences: the flow slowly evolves to a different purely 
hydrodynamical state, and the magnetic field now grows 
exponentially, but both of these effects are sufficiently slow 
for the wavelike character to persist over many wave peri- 
ods. 

Two other interpretations for the nature of these waves 
remain possible: either hydrodynamic fluctuations (e.g. 
inertial waves or Rossby waves) or dynamo waves, as ex- 
pected on the Sun. These possibilities were tested by com- 
paring oscillations in the velocity field and in the magnetic 
field in the kinematic simulations. We found that a high 
frequency signal is present both in the flow and in the 
magnetic field. This demonstrates the presence of hydro- 
dynamic waves, which induce magnetic fluctuations. The 
lower frequency signal is however absent in the flow. This 
provides a proof of their "dynamo wave" nature. 

We have numerically observed such dynamo waves for 
aspects ratio up to 0.8. For the parameters investigated 
here, the transition from a dynamo dominated by a 
fluctuating dipole to a dynamo wave occurs for an aspect 
ratio close to 0.65. This transition exhibits hysteresis: 
once a dynamo-wave solution is present, the aspect ratio 
can be reduced again down to 0.6, while maintaining this 
dynamo mode. 
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Fig. 2: Radial magnetic field at the surface of the outer sphere, 
for aspect ratios of 0.45 (a) and 0.65 (b). Azimuthal magnetic 
field below the surface of the 0.65 aspect ratio model (c). 



Connections with parameterized models. But- 
terfly diagrams indicative of the solar cycle are usually 
produced using simplified parameterized models or "mean 
field" models. These models require a prescription of the 
turbulent induction, the so-called "a-effect" (which can 
also be introduced in terms of deviation from axisymme- 
try [12]). We should stress that this is a valid approxi- 
mation only if certain conditions are satisfied (e.g. [13]). 
Such butterfly-like diagrams are generally not produced 
by direct three-dimensional modelling, with the notable 
exception (only in the reverse direction) of the pioneering 
work of Gilman and Glatzmaier [5,14]. 

Because of the strong symmetry of the convective flows 
influenced by the rapid rotation of the planet or the star, 
it is well known that two independent families of solu- 
tions exist, namely with dipole symmetry (antisymmetric 
with respect to the equator) and quadrupole symmetry 
(symmetric with respect to the equator). Both families 
of solutions are often described in reduced parameterised 
models [15, 16] , and we have observed these two families in 
our fully 3D simulations (Fig. 3). Both branches are stable 
in our simulations for long periods of time, but can also be 
destabilised to yield a change of symmetry. In fact, despite 
the relative complexity of our model, the temporal behav- 
ior of both symmetries is clearly reminiscent of kinematic 
studies of earlier reduced models (Fig. 4c, d and [15]). 

The simpler meanfield equations for the axisymmetric 
field are obtained by writing the flow and field as 



u = sujeri 



B = B p + Be = Vx(Ae )+£e , (6) 



i.e. assuming a mean flow in the form of a zonal shear 
only. In the isotropic case, the axisymmetric part of (2) 
yields (e.g. [3]) 



dA 

— =aB + Rm- 1 V 2 A, 
at 



(7) 



dB 

— = s B„ • Vw + (V x aBJ • + Rm" 1 D 2 B , (8) 
at 

where s denotes the cylindrical radius and T>2 = A — 1/s 2 
(note that Pm in (2) is here changed to Rm as the flow is 
now assumed to be given). 
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Fig. 3: Time evolution of the zonal average of the azimuthal 
magnetic field below the surface of the model, for an aspect ra- 
tio of 0.65: the antisymmetric (a) and symmetric (b) solutions. 



For an instability of (7-8) to exist, these equations must 
not decouple (this is the essence of Cowling's anti-dynamo 
theorem [17]). Equation (8) involves A through two terms. 
Reduced models have been classified in two categories de- 
pending on the dominant term. The first term on the RHS 
of (8) involves the zonal shear and is referred to as the f2- 
effect. The second term in the RHS of (8), as well as the 
first term on the RHS of (7), involve mean induction from 
non-axisymmetric features in the flow and are referred to 
as the a-effect. 

Dropping the a-effect term in (8) and writing the re- 
sulting equations in a simplified cartesian geometry yields 



dA _ _^ i . . dB 

— =aB + Rm AA , — 
at at 



Ft A 

G^+Rm^AB, (9) 
ox 



where G = du y /dz. Parker [18] was the first to identify 
travelling waves solutions (dynamo waves) of the above 
system. These oscillatory dynamos, named Parker waves, 
were obtained by Roberts [15] for nearly axisymmetric dy- 
namos in spherical geometries (following the formalism of 
Braginsky [12]). It was found that while the afi-dynamos 
tended to be oscillatory (complex growth rate), for a 2 - 
dynamos the simplest dipole solutions tended to be sta- 
tionary (real growth rate). A similar behavior can eas- 
ily be traced in the simpler cartesian example above (see 
also [19,20] for a discussion of the generic behavior of such 
nearly axisymmetric mean field dynamos). 

We can perform further comparisons with reduced 
models by studying only the axisymmetric component 
of the simulated field. Figure 4 shows the azimuthally 
averaged field for some of our fully 3D simulations. 
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Fig. 4: The zonal average of the magnetic field in our 3D simulations. Contours of the toroidal (east-west) part of the field are 
plotted in the left hemisphere and lines of force of the meridional (poloidal) part of the field plotted in the right hemisphere. 
The aspect ratio is increased from 0.45 (a) to 0.6 (b) and to 0.65 (c-d). The sequence of dynamo waves is represented for the 
antisymmetric mode (c) and symmetric mode (d). It is similar in nature to that produced by parameterized models [15]. 
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The Earth-like mode is represented for aspect ratios 
of 0.45 and 0.6 (a & b). The active dynamo region 
lies outside the tangent cylinder [4], it therefore gets 
increasingly constrained as the inner-sphere in increased. 
The dipole eventually drops for large aspect ratio, when 
the volume outside the tangent cylinder becomes too 
small. Weakly dipolar solutions were also obtained at 
large aspect ratio in simulations using equations modified 
by hyperviscosity [21]. The dipolar solution was also 
found to decay and eventually vanish by increasing the 
aspect ratio in a reduced parameterized model for the 
Earth's core [22]. Here we show that the steady dynamo 
branch can be replaced, at larger aspect ratio, by an 
oscillatory dynamo mode. Comparison with reduced 
parameterized models can help interpret this transition 
to the solar-like mode. A strong zonal wind develops, in 
our simulations, in the Solar-like mode. Although the 
terminology of parameterized models must be used with 
care for direct simulations (the hypothesis of scales sepa- 
ration does not strictly apply), this suggests a transition 
from a dynamo of the a 2 type to a dynamo of the a£l 
type as the aspect ratio is increased. Indeed, Earth-like 
three-dimensional models have been interpreted in terms 
of regeneration by convective vortices only, and thus 
closer to the a 2 formalism [23] (sometimes referred to 
as "giant a-effect"), whereas the a£l formalism provides 
the classical framework to model solar dynamo-waves, as 
guided by the strong shear at the base of the convection 
zone [15, 16]. Such nearly axisymmctric dynamos [15] 
produce cyclic magnetic behaviour very similar to the 
cycles examplified on Figure 3. 

Conclusions. - By varying the aspect ratio, we have ob- 
served a sharp transition from a dipole dominated large 
scale-magnetic field to a cyclic dynamo with a weaker 
dipole. This indicates that the geometry of the dynamo 
region sevcrly constrains the existence of the dipole dom- 
inated solution. We should however stress that other pa- 
rameters, involving ratio of typical forces, could affect the 
precise value of the critical aspect ratio for transition. The 
values of these parameters in our simulations (as in all nu- 
merical models to date) are indeed very remote from the 
actual relevant values for the Sun, or for the Earth. The 
potentially strong effect of this parameter change on the 
dynamo solution should not be under-estimated. It is in- 
deed quite striking, that despite these shortcomings, nu- 
merical models can capture a good part of the qualitative 
fcathures of the solar and geo-magnctic fields. 

Recent observations of stellar magnetism appear to cor- 
roborate this mechanism. Donati et al. [1] reported ob- 
servations of a strongly dipolar field in a fully convective 
star (V374 Peg). More recently, Donati et al. [24] report 
magnetic observations of r-Bootis, a rapidly rotating F 
star, i.e. one with a relatively shallow outer convection 
zone. Not only did they observe a rather complex mag- 
netic field structure, but they also report that the overall 
polarity of the magnetic field has reversed after one year 



of observation. They interpreted this observation as an 
indication that the large aspect ratio r-Bootis star is un- 
dergoing magnetic cycles, similar to those of the Sun. 

Futhcr observations of planets and stars arc needed, but 
clearly the observations available so far seem to confirm 
the important role of the aspect ratio in controling the 
transition from steady to cyclic dynamo modes. 

* * * 
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